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Abstract 

Motivated by recent results from SuperKamiokande, we study both solar and atmo- 
spheric neutrino fluxes in the context of oscillations of the three known neutrinos. We 
aim at a global view which identifies the various possibilities, rather than attempting 
the most accurate determination of the parameters of each scenario. For solar neu- 
trinos we emphasise the importance of performing a general analysis, independent 
of any particular solar model and we consider the possibility that any one of the 
techniques — chlorine, gallium or water Cerenkov — has a large unknown systematic 
error, so that its results should be discarded. The atmospheric neutrino anomaly is 
studied by paying special attention to the ratios of upward and downward going v e 



and Up fluxes. Both anomalies can be described in a minimal scheme where the re- 
spective oscillation frequencies are widely separated or in non-minimal schemes with 
two comparable oscillation frequencies. 

We discuss explicit forms of neutrino mass matrices in which both atmospheric and 
solar neutrino fluxes are explained. In the minimal scheme we identify only two 
'zeroth order' textures that can result from unbroken symmetries. Finally we discuss 
experimental strategies for the determination of the various oscillation parameters. 



^This work was supported in part by the U.S. Department of Energy under Contracts DE-AC03-76SF00098, in part by the 
National Science Foundation under grant PHY-95-14797, in part by the TMR network under the EEC contract n. ERBFMRX- 
CT960090. 



1 Introduction 



The solar and atmospheric neutrino flux anomalies have both been considerably strengthened by recent observa- 
tions from Super-Kamiokande. The solar neutrino flux is measured to be jjj 0.37±0.03 of that expected from the 
'BP95' standard solar model ||, without including any theoretical error. This is the fifth solar neutrino exper- 
iment to report results in strong disagreement with the predictions of solar models. Furthermore, using a solar 
model independent analysis, the measured solar fluxes are found to be in conflict with each other. For events at 
SupcrKamiokandc with visible energies of order a GeV, the ratio of 1 ring ^t-like to e-like events is 0.66 ± 0.10 
that expected from calculations of the flux of neutrinos produced in the atmosphere in cosmic ray showers || . 
Furthermore, the distribution in zenith angle of these 1 ring events provides striking evidence for a depiction of 

which depends on the distance travelled by the neutrinos before reaching the Super-Kamiokande detector. 
In particular, the observed up/down ratio of the multi-GeV, /x-like events is 0.52 ± 0.07. This significantly 
strengthens the evidence that oscillate as they traverse the earth. 

In this paper, we interpret the solar and atmospheric neutrino flux anomalies in terms of oscillations of the 
three known neutrinos v ellT . The lightness of these three neutrinos, relative to the charged fcrmions, can be 
simply understood as resulting from large SU(2)^ ® U(l)y invariant masses for the right-handed neutrinos, via 
the see-saw mechanism. We do not consider the possibility of a fourth light neutrino, as it would have to be 
singlet under SU(2)i ® U(l)y, and would either require a new mass scale far below the weak scale, running 
counter to the idea of the see-saw mechanism, or a more complicated see-saw. 

Theoretical ideas about generation mixing are guided by the quark sector, where the mixing angles are all 
small, indicating a hierarchical breaking of horizontal symmetries in nature. A similar hierarchy of horizontal 
symmetry breaking in the lepton sector is also likely to yield small angles, suggesting small probabilities for a 
neutrino to oscillate from one flavour to another. However, the solar and atmospheric neutrino flux measure- 
ments both require neutrino survival probabilities, P ee and -P w , far from unity. Over a decade ago [Q, it was 
realised that large angles were not necessary to account for the large suppression of solar neutrino fluxes - 
while v e have charged current interactions in the solar medium, za, T do not, allowing a level crossing phenomena 
where a u e state produced in the solar interior evolves to a v^ T state as it traverses the sun. This simple picture 
can reconcile the three types of solar neutrino flux measurements with the standard solar model, for a mixing 
angle as small as 0.03 — a significant achievement. Could such resonant oscillations occur for atmospheric 
neutrinos in the earth, again allowing a small vacuum mixing angle? In this case, since the earth docs not have 
a continuously varying density, the matter mixing angle in the earth is much larger than the vacuum mixing 
angle only in a small range of energies. Hence, an oscillation interpretation of the atmospheric neutrino fluxes 
requires a large mixing angle, and calls into question the frequently stated theoretical prejudice in favour of 
small mixing angles. 

In this paper, we attempt to understand both solar and atmospheric neutrino fluxes using 3-generation 
neutrino oscillations, aiming at a global view which identifies the various possibilities, rather than attempting 
the most accurate determination of the parameters of each scenario. When data from chlorine, gallium and 
water Cerenkov detectors are fitted to a standard solar model, standard analyses find very small regions of 
neutrino mass and mixing parameters. For 2-gcneration mixing, these arc known as the "small angle MSW" , 
"large angle MSW" and "just so" regions. This analysis has been extended to the case of three generations ||, 
with a single matter resonance in the sun, as suggested by the atmospheric neutrino data. The large and small 
angle MSW areas are found to merge into a single MSW volume of parameter space. In subsection 2.1, we 



study how this volume is enlarged when a solar model independent analysis of the solar fluxes replaces the use 



of a single solar model. In subsection 2.2 we extend our analysis to see what areas of neutrino parameter space 
become allowed if one of the three observational techniques to measure the solar fluxes is seriously in error. 

We combine these regions of parameters with those yielding the atmospheric fluxes, and find there is still 
considerable allowed ranges of masses and mixing angles. This is done in section |3|, assuming that the smallest of 
the two neutrino squared mass differences is too small to affect the oscillations of atmospheric neutrinos (minimal 
scheme) . In section ^, on the contrary, we allow for the possibility that the two independent neutrino squared 
mass differences are both large enough to affect atmospheric neutrino oscillations (non minimal schemes). For 
solar neutrinos, this requires that there is a serious flaw either in at least one measurement technique or in solar 
model analyses. 

The forms of neutrino mass matrices that can lead to a large ^ v T mixing for atmospheric neutrinos 
are discussed in section |5[ In section |(] only two 'zeroth order' textures for neutrinos masses are identified that 
can account for the atmospheric and solar neutrino data in the minimal scheme and can result from unbroken 
symmetries. 

Our conclusions are drawn in section |^. Based on a simple set of alternative hypotheses, we discuss how 
future measurements could eventually determine the two neutrino mass differences and the three mixing angles. 
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2 Solar neutrinos: model-independent analysis 



In the flavour eigenstatc basis, in which the charged lepton mass matrix is diagonal, the neutrino mass matrix 
is in general non-diagonal. It may be diagonalized by a unitary transformation: 

VJ = VfiUt (2.1) 

where Vf and vi are flavour and mass eigenstate fields, respectively. The leptonic analogue of the Cabibbo- 
Kobayashi-Maskawa mixing matrix is V T , since the W boson couples to the charged current v^V^^e^. In 
addition to the three Eulcr angles, V contains physical phases: one if the light neutrinos are Dirac, and three 
if they are Majorana. These flavour and mass eigenstate fields destroy basis states which are related by 

Wf)=V f M) (2.2) 

If some process creates a flavour eigenstate, \vf), at time t = 0, then at a later time t it will have evolved to the 
state \vf,t) = ipf(t)\i>fi} via the matrix Schroedingcr equation 

i^ = {V^V* + Acc + £)1> (2.3) 

where E is the energy of the relativistic neutrino, m„ is the diagonal neutrino mass matrix with entries m,, £ is 
an irrelevant term proportional to the unit matrix, and Aqq represents matter effects. For neutrinos propagating 
in matter with electron number density N e , Aqq is a matrix with a single non-zero entry, Aq C = \f2GpN e . 
The mixing matrix V can be written quite generally as 

0\ /l 

V = i? 23 (M | e** Ru(9 13 )R 12 (e 12 ) e ia (2.4) 

I J \0 e if} 

where Rij(9ij) represents a rotation by 6{j in the ij plane. We have chosen a sequence of rotations which 
frequently arises in the diagonalization of simple hierarchical forms for the neutrino mass matrix, as illustrated 




in section |6J. From equation (2.3) we see that the phases a and (3 never appear in oscillation phenomena, and 
hence can be dropped, giving 

C12C13 C13S12 
V = ( -c 23 si 2 e^ - C12S13S23 ci 2 C23e^ - S12S13S23 ci 3 s 23 | . (2.5) 
S23Sl2e 10 - C12C23S13 -ci 2 S23e 10 - C23S12S13 

Each Rij must diagonalize a symmetric 2x2 sub-matrix determining tan 2(9^ , hence, without loss of generality, 
we may choose < Oij < 7r/2, while < <f> < 2tt. A more convenient choice is to keep #12, 13 in the first quadrant, 
while < #23,</> < 7T- We choose to order the neutrino mass eigenstates so that Ato^ > Amf 2 > 0, where 
Amy = mf — rrij. Notice that with this parametrization V e 3 -C 1 means ^13 close to or to 90°. 

To study solar neutrinos, we are interested only in the electron neutrino survival probability, P ee , and hence 




in the evolution of tp e . This evolution does not depend on #23 or on — on substituting (2.4) in (2.3), -R23 and 
4> can be absorbed into redefined states [tl and r' . Hence, we have shown quite generally that P ee depends only 
on four neutrino parameters: Amf 2) Am| 3 , #12 and #13. 

For an oscillation explanation of the atmospheric neutrino fluxes, Am^ is sufficiently large that it does not 



cause a resonance transition in the sun. In the Landau-Zehner approximation, the evolution equation (2.3) can 
be solved to give M 



(i-p p o\ /iK?r 

Pee = (|Kl| 2 , |K2| 2 , |K 3 | 2 ) P 1 - P |K2 I 2 I (2.6) 

V 01/ VlKTl 2 

where V™ are the mixing matrix elements in matter, and P is the transition probability between the states at 
resonance: 

P_„-£Wis fl /p F\ F - ^ ATO i2sin 2 (2gi 2 ) Aw i2 cos20 12 

0{E - Ea) ' ^"41^^1,003(2^)' EA -2V2G F \N e \oco S ^ 13 ^ 

Here E is the neutrino energy, 9 is the step function, the 1 subscript indicates that N e and its gradient dN e /dx 
are evaluated at the resonance point, while the subscript indicates the production point. The large mass 
splitting ATO23 enters P ee only via the matter mixing angles, and decouples from these expressions in the limit 
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Flux of boron neutrinos, $sb, in 10 6 / cm 2 sec 

Figure 1: Values o/($8 B ,$tb ) measured by the Chlorine experiment (continuous lines), the Gallium experiment 
(dashed lines) and by the SuperKamiokande experiment (long dashed lines) assuming various neutrino oscillation 
schemes: • no oscillation in fig. |Ja; • an energy-independent P(y e — > u e ) = 0.85 in fig. ^p: • the best-fit point 
of the small-angle MSW oscillation in fig. Hie; • the best-fit point of the large-angle MSW oscillation in fig. Ihi. 



that it is much larger than AnE, and also in the limit that #13 vanishes. For most of this section we make 
Am^ sufficiently large that it decouples, and we comment at the end on the effect on the allowed regions of 
parameter space for non-zero #13 and small Am 2 ^, where Am^ effects may not decouple. 
The signals Si at the three types of solar neutrino experiments are 

Si = JdE^(E)[at(E)P ee {E)+af(E)(l-P ee (E))}, z = {SK,Ga,Cl} (2.8) 

where &(E) is the total flux of solar neutrinos with energy E, and a1^(E) are the interaction cross sections at 
experiment i for electron-type and non-electron-typc neutrinos, respectively (only the water Cerenkov detectors 
are sensitive to neutral currents, so a & a {E) = <j^{E) = 0). We will use the theoretical predictions of the various 
cross sections found in @| . The flux $(£) is broken into components in the standard way by specifying the 
production reaction, giving j?J 

/>oo 

= X>a/a(£), With / f a (E)dE = l (2.9) 

and a = pp, pep, 7 Be, 13 N, 15 0, 17 F, ^ hep. At this point we follow the (nearly) model-independent treatment of 
the fluxes described in || by making the following assumptions: 

1. The energy dependence f a {E) of the single components of the neutrino fluxes predicted by solar models 
(H, |^| for instance) are correct. In fact the f a {E) do not depend on the structure of the sun, and are the 
same in any solar model that does not introduce non-standard electroweak effects H . 

2. The overall <f> Q can differ from their solar models predictions. However there are strong physical reasons 
to believe that the ratios <3?i3n/<I>i5o and <f> pe p/<5pp can be set to their solar SM values Furthermore we 
neglect entirely hep and 17 F neutrinos, which we expect to be extremely rare. 

3. The present total luminosity of the sun, Kq, determines its present total neutrino luminosity as 

ck ^ a 

where Q = 26.73 MeV is the energy released in the reaction 4p + 2e — > 4 Hc + 2is e , and Kq = 8.53 ■ 10 11 



MeV cm -2 s _1 is the solar radiative flux at the earth. Using ( 2.1C ) amounts to assuming that the solar 
energy comes from nuclear reactions that reach completion, and that the sun is essentially static over the 
10 4 years employed by photons to random-walk out of the solar interior. 
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Figure 2: Values of the x-sguare 
as function of an energy independent 
P{y e — ► h> e ). The parameter A is de- 
fined in eq. ( [2.141 ). Also shown is the 
X 2 with one experiment discarded and 
A = 1. 



After the first assumption we have one free parameter $ a for each a; the second then reduces the number of 
free parameters to four, which we can take to be 



<I>, 



$pp + $pep, 



$CNO = $13N + 



r>0 . 



T3c 



and 



$8 



(2.11) 



The luminosity constraint allows us to eliminate <& p , giving 



Si = 5i(Amf 2 , 6*12, 6»i 3 ;$8 B ,$7B e 



^CNO 



(2.12) 



Since solar models give a stable prediction for <E'cno/ < £' 8 b = 0-22 we have singled out this ratio and we will 
use its SSM value in our analysis. Variations of even an order of magnitude in the ratio affect negligibly our 
final results, since the two neutrino components have similar cross sections in existing detectors. 



2.1 Model-independent solar analysis — all experiments 

The signals now depend only on <1>sb and •J'tbcj so that, for any given oscillation pattern P ee (E) it is possible to 
plot the three experimental results^] 



ocxp 

°C1 — 

oexp 

QCXp 



-36 s -l 



(2.54 ±0.20) 10" 
(75 ±7) 10" 
(2.51 ±0.16) • 10 B cnT's 



(2.13a) 
(2.13b) 
(2.13c) 



as three bands in the (^sbj $tb c ) plane. The three bands will in general not meet, giving interesting solar model 
independent restrictions on the oscillations parameters. 

We begin the analysis by studying the case of no neutrino oscillations (P ee = 1). In this particular case 
the solar model independent analysis does not give a strong result. Surprisingly the three bands perfectly 
meet [l3| as shown in fig. [l]a, but mainly in the unphysical &TB e < region, with a small area in the physical 
region lying within 2a of each central value. Since the physical crossing region has a negligible 7 Bc flux, the 
value of $cNo/ < I )7 Be becomes completely irrelevant. 

To discuss this case in a quantitative way and to deal with more general cases it is useful to introduce the 
X-square function 



Xl (Pee(Am 2 12 , 6 V 2, 013), $*B, $?Bc) = 



Si - sr 

A£ cxp 



E 



($, 



^SSM 



)(** 



A 2 A$ 2SSM 



(2.14) 



where AS°* P is the la uncertainty for experiment i, given in ( 2.13 ), $ SSM is the flux prediction of the solar 
model H and A<I> SSM is the corresponding error matrix, taken with some generosity. The la ranges of $»b and 

lr The SuperKamiokandc experimentalists give directly the value of the flux they measure. The other experiments involve more 
uncertain neutrino cross sections and prefer to give the frequency of events measured per target atom in their detector. For 
simplicity we have omitted this detail in the text, leaving a trivial inconsistency between eq. (2.13c) and (2. 8). 
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Figure 3: Allowed regions in the plane (sin 2 2#i2, Amf 2 ) for #13 = 0, 15° and 30°. The upper plots assume that 
the BP solar model is correct. The lower plots are the result of the solar model independent analysis described 
in the text. 



<i>7B are represented by the ellipse in fig. [|. We perform our analysis with two choices for A$ = A • A$ SSM . 
We call the first choice, A$ = A$ SSM , ll solar SM inspired"). The second choice, A$ = 8 • A$ SSM ("model 
independent" ) has the same shape as the first, but is eight times as large. The part of the analysis done using 
this A$ is virtually free of solar physics input. The choice A = 8 (rather than A = oo) avoids unnatural 
values of $tbc ■ This choice essentially ignores solar physics considerations, but the virtue of having a number 
of independent experimental results is precisely that we need no longer rely heavily on solar modelling to gain 
insight into the underlying particle physics. 

Minimizing the \ 2 m the positive flux region we obtain minx|(P ee = 1) = 8.25. The usual criterion 
for goodness of fit says that a \ 2 with one degree of freedom larger than 8.25 is obtained with a very small 
probability, p « 0.4% (a careful Monte Carlo treatment of the 3>tbc > constraint gives similar results |l5|]). 
We however remark that, if the sun really emits the best-fit fluxes, $?b c = and $sb = 2.5 10 6 /cm 2 s, there is 
a 10% probability that statistical fluctuations produce the present experimental data. 

We can just as easily investigate the slightly more general case of an energy independent P ee . The dependence 
on the neutrino parameters Am 2 2 , # 12 , and $13 arises through P ee ; if the survival probability is a constant, then 
we can minimize X\ m the positive-flux region for any value of P ee to obtain minx 2 (P ee ), which is plotted in 
fig. U for A = 1 (SSM analysis), A = 8 (SSM independent analysis) and A = 00 (completely model independent 
analysis). For P ee ~ 0.85, minx! drops to 5, but the fluxes required to get relatively small x 2 values are 
disfavoured by solar physics considerations — 'I'cno and ^^Bo must be nearly made to vanish, as shown in 
fig |l|b. When P ee < 1/2 the (accidental?) threefold crossing no longer occurs, so that this case can be fi rmly 
excluded in a solar- model independent way || Q (see fig. |f). However, as we shall see in subsection 2/2, once 



we allow for the possibility that one type of experiment's results should be discarded, it is possible to obtain 
good fits of the data for constant P ee ~ 1/2 without having to resort to unnatural flux values. 
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Figure 4: Isoplot of Xs> minimized in 
the mixing parameters. 



Of course, we are interested in any points in parameter space that fit the data well, regardless of whether 
they lead to constant P ee . For any values of Am? 2 , #12, and 613 we can make plots similar to fig. |l|a. Fig.s 0c 
and [j]d show two examples that illustrate the familiar 2-generation small and large angle MSW solutions, which 
evidently fit the data well if standard solar model fluxes are used. 

In fig. [| we show how the allowed regions in neutrino parameter space change if we let the fluxes vary over 
an expanded range of values. For each point in (A?n 2 2 , 6*12, #13) space, we minimize Xi an d xi by varying 
the fluxes within the physical region, and then we plot contours of minx 2 in the (sin (26*12), Am 2 2 ) plane for 
various values of 6*13. The results for the "SSM inspired" and "model independent" analyses are shown in fig.s || 
(upper row and lower row, respectively). The contours are for \ 2 — 3 and x 2 — 6. 

For small 6*13 the "SSM inspired" results show the standard small and large angle MSW regions. For larger 
values of 613, the two MSW regions join, and, as 6*13 approaches 7r/4, the solutions with large 6*i 2 disappear. 
For 813 = 7r/4 the region with minx? < 3 is in fact absent entirely. 

The "model independent" results similarly exhibit a very strong 6*13 dependence. We see that the "model 
independent" analysis continues to give strong restrictions of the oscillation parameters — in particular the 
Am 2 2 values with minx§ < 3 are always in the range ~ lO^ 4 ^ 5 -* cV 2 . This will not remain true when we 



consider the consequences of ignoring one experiment's data in subsection 2.2. 

If < I ) cNo/ < I )s B is ten times larger than in SSM there are new allowed regions. However these possible new 
regions, with Am 2 = lCT( 5 - 6 > eV 2 and sin 26*12 <^ 10 , are excluded in a model-independent way by the non 
observation of a day/night asymmetry at SuperKamiokandc pi ^5| . The recent data jlj on this asymmetry 
in fact disfavour as well the large angle MSW solution of the SSM-inspired analysis. Moreover, we have not 
included in our x 2 analysis the SupcrKamiokande measurement of the distortion of the ^ spectrum (l], [f5| , 
because the present positive la signal could be produced by a ^hep/'&SB ratio 15 times larger than the prediction 
of BP95 @. Without a very large hep flux, the present measurement excludes an otherwise allowed region with 
Am 2 « 10" 4 eV 2 and sin 2 26> X2 in the range 1CT 4 -j- 1CT 1 111 [HI. 



Our model independent analysis allows us to investigate how well present experiments are able to measure 
the SSM-independent neutrino fluxes $sb and $7Bc- This question is answered in fig. ||, where we plot the values 
of the fluxes that can give a good (x| < 6) or very good (xt < 3) fit for some value of the oscillation parameters 
Am 2 2 , 6*12 and 6*13. We see that the value of <I>8b is currently determined with an error larger than the solar 
model expectation. It will be directly measured in the new on-going SNO experiment. On the contrary the 
value of $7B e is at present totally unknown: in fact in the small angle MSW solution the monochromatic 'Be 
flux can be completely converted into ^ neutrinos, that are not detected by existing experiments. Borexino will 
be able to detect neutral currents effects in this range of energies and probably allow a direct determination of 

As discussed above we perform our analysis under the assumption that Am 23 is large enough that its effects 
decouple. For any given Ato 23 it is straightforward to reproduce fig. [?] by using the exact expressions for 



and 6*5^ in equation (2.6). In this way we find that for small 6*13 ( < 15°), our results are insensitive to Am 23 
down to Am 23 = 5 • 10 _4 eV 2 . For large 6*13, Am| 3 effects start to become noticeable when Am 23 drops below 
2 x 10~ 3 eV 2 ; for example, for 6*13 = 40° and Ato 23 = 5 • 10~ 4 eV 2 , the allowed region in the SM inspired analysis 
is significantly smaller than in the decoupled limit, with the Xmin < ^ region never reaching sin 2 (26 l i2) > 0.1. 
In spite of these changes for small Am 23 , the essential features of fig. |3| in any case remain unchanged. 
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Figure 5: fits of the solar data in the plane (sin 2 2#i2, Am 2 2 / cV 2 ) for #13 = 0, 15° and 30° assuming that one of 
the three solar neutrino experiments has a large unknown systematic error ( Super Kamiokande in the first row, 
Chlorine in the second and Gallium in the third) and is therefore discarded from the analysis. The contours are 
for x 2 = 3 and \ 2 = 6. 
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2.2 Model independent solar analysis — one experiment ignored 

In subsection ^T|, the present level of experimental evidence allowed us to omit one restriction (the solar model) 
and still yield interesting results. Likewise, we can choose to omit one experiment from the analysis while 
keeping some solar information and still yield interesting results. 

The motivation for this is obvious: neutrino experiments are extremely difficult to perform and particular 
detection schemes may suffer from some systematic error previously not considered. We make no judgements 
here about the errors associated with any particular experiment. Instead we consider analyses where we do not 
include one class of experiment, either water-Cerenkov, gallium or chlorine, which we designate 0i and 0a 
respectively. However, because we are losing an experiment, it is impossible to analyze the data without some 
level of information regarding the solar model. Consequently, we perform the analysis only within the solar 
SM inspired region. The results for this analysis are shown in figures |5| (upper row: without SuperKamiokande 
data, middle row: without chlorine data, and lower row: without gallium data). 

The S'ff case largely resembles the complete data set analysis, with some additional space allowed in the 
higher Am 2 2 region. In contrast, the other two cases (0a and 0V) show considerable differences. 

For the 0a case, there is a strong preference for either small Ato 2 2 or large 6*12 and 6*13. For the 01 case, 
for both large #12 and 6*13 we have the presence of large regions with large Am 2 2 = 10~ 4 eV 2 , above the level- 
crossing threshold, and with small Ato 2 2 , m t ne non-adiabatic region. In either case, in a large portion of these 
regions matter enhancements are unimportant. That is, in the absence of one of these two classes of experiment, 
given sufficiently large angles, the solar neutrino problem can be resolved simply by vacuum oscillations alone! 
In such a case, new experiments, such as Borexino, would see an absence of energy dependence in the electron 
neutrino survival probability. 



3 Atmospheric and Solar Neutrinos: The Minimal Scheme 

The simplest picture for reconciling both solar and atmospheric neutrino fluxes via oscillations of v e .ii,T results 
when there is a hierarchy |Am 23 | 3> |Am 2 2 |, and Ato 2 2 is too small to affect oscillations of atmospheric 
neutrinos. In section ||, we showed that in this case the solar fluxes depend only on Ato 2 2 , 612 and #130, and 
below we show that the atmospheric fluxes depend only on Am 23 , 9 23 and #13. In the limit that #13 = 0, the two 
phenomena become independent, in the sense that they depend on no common parameters: solar oscillations are 
v e — > Vp, at a low frequency, while atmospheric oscillations are — > v T at a much higher frequency. However, 
solar oscillations are allowed for a wide range of parameters with large #13 , and the atmospheric data does not 
require #13 to be very small. Hence, in this section we explore this simple picture keeping #13 as a free parameter. 
We comment on the alternative possibility — that Ato 2 2 is large enough to contribute to atmospheric neutrino 
oscillations — in section [|. 

Matter effects in the earth are important only for a relatively sm all f raction of the atmospheric neutrinos, 
those with high energy, and they are neglected heref]. In this case, (|2.3| ) can be integrated to give oscillation 
probabilities Pff(t) = \Affi(t)\ 2 , where A is given by the matrix equation 

A(t) = Ve- iEt V^. (3.1) 

Since an overall phase in A is irrelevant to P, and Ato 2 2 effects are negligible, we may make the substitution 

e -iEt diag ^ 1; e -iAmL*/2B) (3 2 ) 



Using the form (2.4) for V, we immediately discover that the probabilities are independent of 6*12 and 0, as well 



as a and (3. The probabilities are given by 

P eM = s 2 3 sin 2 20 13 S 23 (3.3a) 

P eT = c\ 3 sin 2 26 13 S 23 (3.3b) 

P\it = 4 sin 2 2023 ^23 (3.3c) 



or equivalcntly, by unitarity 



P ee = 1 - sin^ 29 13 S 23 (3.3d) 
= l-4c 2 3 .s 2 3 (l- c 2 13 sl 3 )S 23 (3.3c) 

Prr = 1 - 4C 2 3 C 2 3 (1 - C 2 3 C 2 3 ) S 23 (3.3 f) 



2 Although for non-zero 6*13, there is a dependence on Am^ if it is small enough. 
3 For more details see e.g. ref. JUT. 
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where = sin 2 (Am| 3 t / '4E) . The parameter Am^ can be extracted from the data by fitting to the zenith 
angle distribution of the events. Here we concentrate on the determination of the parameters #13 and 823- These 
can be extracted, independent of the value of Ato 2 3 , if we assume that the downward going neutrinos have not 
oscillated, while the upward going neutrinos are completely oscillated, so that 6*23 is averaged to 0.5. In view of 
the reported angular distribution of the multi-GeV data for 1-ring e-like, 1-ring p-like and partially contained 
(PC) events ||, this assumption appears to be valid, at least for angular cone sizes about the vertical which are 
not too large. For events of class i, which are induced by v e charged current, za, charged current and neutral 
current interactions with relative probabilities /* cc , P^cc anc ^ /Arc ^ ne up-down ratio pi is given by 

7V T 1 
ft = = flee ■ (Pee + rP e „) + flee ■ (P^ + -P^) + Pnc- (3-4) 

T I 

where we have set 5*23 = 0.5, and N± are the number of upward and downward events of class i. We are 
interested in i being 1-ring e-like, 1-ring /x-like and PC. The overall normalization of these event numbers has 
considerable uncertainties due to the calculation of the neutrino fluxes produced in cosmic ray showers, hence 
we consider three up-down ratios 

p e = 1.23 ±0.29 (3.5a) 
= 0.62 ±0.16 (3.5b) 
pre = 0.48 ±0.12 (3.5c) 



and two ratios of downward going fluxes 



N, 



£r = 3.0 ±0.6 (3.5d) 



— = 1.3 ±0.3 (3.5e) 

where r is the ratio of to v e fluxes. The numbers give the Super-Kamiokande data, extracted from the figures 
of Rcf. jj), with upward and downward directions defined by the azimuthal angle having cos8 within 0.4 of the 



vertical direction. The parameters £r and £' represent the theoretical values for the ratios of ( 3.5d) and (3.5e). 
These two downward going ratios do not involve oscillations, and the Super-Kamiokande collaboration compute 
Monte Carlo values of 3.1 and 1.0, respectively, agreeing very well with the data. Since these two ratios do 
not probe oscillations, at least within our assumptions, we do not use them for the fits below. We do not use 
the sub-GeV data as the poor angular correlation between the neutrino and charged lepton directions leads 
to a smoothing of the up-down ratio. From the flux calculations of Honda et al |18| , and using the measured 
momentum distributions for the events ||, we estimate r = 4.0 ± 0.5, for this multi-GeV data near the vertical 
direction. A more refined analysis would use a larger value of r for PC events than for FC events. 

The results of a fit of the three up/down ratios to the two free parameters 023 and #13 are shown in 
figure ||(a). We have obtained the fractions P eCC mCC nc from the Mon te Ca rlo results of the Super-Kamiokande 



collaboration |]3j, and we have used the oscillation probabilities of ( 3.3c ). In order to work with Gaussian 

T I 

distributed experimental data, we have directly fitted the six measured neutrino numbers N- leaving arbitrary 
the three overall fluxes of each type, Nj + N^. The preferred region of the plot is easy to understand, since at 
the point 823 = 45° and # 13 = 0, the u e are unmixed, while there is complete i/„ <-> v T mixing, so p e ~ 1 and 
p^ ~ ppc — 0.5. It is apparent from Fig. ||(a) that this minimal scheme is allowed for a large range of angles 
about this point: 8 23 = 45° ± 15° and 8 13 = + 45°. 

If the solar neutrino fluxes, measured by all three techniques, are to agree with solar model inspired values, 
then the results of section show that Am 2 2 is too small to affect atmospheric oscillations, it is either of order 
10~ 4 - 10~ 5 eV 2 or of order 10~ 10 eV 2 . In this case, the minimal scheme for atmospheric neutrinos, described 
in this section, is the unique possibility using just the three known neutrinos. This observation enhances the 
importance of the fit of figure ^(a); further data will reduce the allowed region, as the three up-down ratios of 



( 3.5c ) have small systematic uncertainties and are statistics limited. The solar neutrino fluxes do not put extra 
constraints on the value of 6*13, although it becomes correlated with 812, as shown in figure ||. If the atmospheric 
flux measurements require Am^ > 2 x 10~ 3 cV 2 , then the limit on P ee from the CHOOZ experiment ||l9| 
requires #13 < 13°. 

Recent analyses [|o| of SuperKamiokande data that make use of MonteCarlo predictions for the angular 
and energy distributions of the atmospheric neutrinos get more stringent constraints on the neutrino oscillation 
parameters. Our fit uses only those data — the ratio of upward and downward multi-GeV neutrinos (the ones 
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in bins 1 and 5 of the angular distribution in ||) — that do not depend on the spectrum of the atmospheric 
neutrinos nor on the precise value of Am 2 , assuming a full averaged oscillation in between. Since statistics gives 
presently the dominant error, we obtain weaker constraints than in |20[ . If we knew that the neutrino mass 
difference relevant for atmospheric neutrinos were close to the center of the presently allowed region, we could 
add to the data to be fitted the intermediate bins 2 and 4 of || (the bins that contain 'oblique' neutrinos). 
We cannot use in any case the multi-GeV data in the intermediate bin 3, that contains 'horizontal' neutrinos. 
Having doubled the statistics, we would find the more stringent contours shown in fig. ^]a', bl', b2'. We remind 
the reader that our fit does not include matter effects 0. 

4 Atmospheric and Solar Neutrinos: Non-Minimal Schemes 

In this section, we study atmospheric neutrinos when two conditions apply. 



The smallest mass splitting is large enough to affect atmospheric neutrino oscillations: Ato 2 2 > 3 x 
10 _4 eV 2 . For solar neutrinos, this implies that there is a serious flaw either in at least one measurement 
technique, or in the solar models. 

The mass splittings are hierarchical Am^ 3> Am 2 2 . This is a simplification, which we relax at the end 
of the section. It includes the interesting possibility that Am^ is large enough to induce the apparent 
oscillations reported by the LSND collaboration pi]] , while Am 2 2 effects are causing both solar and 
atmospheric oscillations. 



Using (2.E), the v e survival probability is 

Pee = 1 - 4 sin 2 20 12 5i2 - s\ 2 sin 2 26 13 S 23 - c 2 2 sin 2 29 13 S 31 (4.1) 
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where Sij = sin 2 (AmfA/AE) . The above two conditions imply that Am^ > 2 x 10 3 eV 2 , so that, for the 
CHOOZ experiment, (gjjj) should be used with S23 = S31 = 0.5. The CHOOZ limit, P ee > 0.9, then gives 
6*i3 < 0.23. If Am 2 2 were also greater than 2 x 10 _3 cV 2 , then for the CHOOZ experiment one also has 
S12 = 0.5, so that 6*12 < 0.23. However, in this case the survival probability for solar neutrinos is the same 
as for the anti-neutrinos at CHOOZ: P ee > 0.9. Hence, given our two conditions, the observed solar neutrino 
fluxes require Am 2 2 < 2 x 10~ 3 eV 2 . 

It is frequently stated that the three known neutrinos cannot explain the LSND, atmospheric and solar 
neutrino anomalies, as this would require three Am 2 with different orders of magnitudes. However, this argument 
no longer applies in the case that either a solar neutrino measurement technique or solar models are incorrect, 
when a single Am 2 could give both atmospheric and solar anomalies. Hence, we consider first the case that 
Am§3 is l arge enough to explain the observations of LSND. The oscillation probabilities induced by Am 2 3 are 
given by ( [3. 3c ). From the limit on P ee from the Bugey reactor, one then concludes 



Am 2 3 > 0.2 eV 2 , (4.2a) 

and 

6»i3 < 0.1 (4.2b) 

which is significantly stronger than the CHOOZ limit. A second possibility, 6*13 close to 90° does not allow any 
significant oscillations of v e and is thus not acceptable to explain the solar neutrino anomaly at a relatively large 
frequency For atmospheric neutrinos, both upward going and downward going, one may then use oscillation 
probabilities with #13 = 00, and with S 23 and S\ 3 both averaged to 0.5: 

P^ = c 2 3 sin 2 20 12 S 12 (4.3a) 

P eT = s 2 3 sin 2 20ia 5*12 (4.3b) 

Ppr = ~\ sin 2 2023 sin 2 20 12 S 12 + \ sin 2 20 23 (4.3c) 

or equivalently, from unitarity 

P ee = 1 - sin 2 2(9 12 S*i2 (4.3d) 

Pw = 1 - \ sin 2 2023 - 4, sin 2 20 12 S 12 (4.3c) 

P TT = 1 - i sin 2 2023 - 4a sin 2 20 12 S 12 . (4.3 f) 

Since in these formulae S\ 3 = 523, 012 enters only via sin 2 20i2 so that, without loss of generality, we may reduce 



uncertainties. We calculate them approximately, using ( |3.4|) with fji cc = f^cc = fuCC = 1 ano - a ^ other 



the range of 9i 2 to < 0i2 < 7r/4. We again study the up- down ratios (3T), as they have small systematic 

fP- _ fPC 
IfiCC ~ Jp.CC 



/-factors equal to zero. A fraction, pf£}} = 1 — sin 2 2023/2, of the downward going oscillate to v T before 



up 

detection, so the up-down ratios are given by 



Pee + rP eil (4.4) 



and 



Hence we find 



MM 

0V - 1) « ~l z ^\ nn ■ (Pe - 1). (4.6) 



r 1 - i sin 2 20 23 



For the multi-GcV data, where the angular correlation is best, r is large, and (4.6) implies that — 1| < 



(l/3)|p e — 1|, in strong disagreement with data of (3.5c). The same inequality holds if p M is replaced by ppci 
when the disagreement with data is even stronger]^ With oscillations of the three known neutrinos, the LSND 
observation conflicts with the atmospheric and solar neutrino anomalies even using a model independent analysis 
of the solar neutrino fluxes or allowing for a systematic error in one of the solar neutrino experiment^. 

Does the atmospheric neutrino data allow other values of Am 23 3> Am 2 2 ? The limit from the Bugey reactor, 
(4.2b), applies for all Am 23 > 0.06 eV 2 , and the up-down ratio relation, (fOj), applies for all Am 23 > 0.1 eV 2 . 



4 In which case the Pij are independent of (f>. 
5 Even ignoring p e , we find ppc,p > 0.61. 

6 For 823 = 0, this corresponds to purely i/ M — » u e oscillations, which is therefore excluded as an explanation of the atmospheric 
neutrino measurements. 
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Hence, A77J23 > 0.1 cV z is excluded. For Am^ < O.leV 2 , the downward going have not oscillated to v T 



when they reach the Super-Kamiokande detector, so that (4.5) is replaced by 



« P m + - r P eil = \-\ sin 2 20 2 3 - ^(p e - 1). (4.7) 



Consistency with the data, (3.5c), is now possible, and requires large #23- As Am| 3 drops below 0.06eV 2 , 
the limit from the Bugey reactor on # 13 is progressively weakened, so that #13 terms must be kept in Py. 
Furthermore, as Am^ drops below 0.01 cV 2 , our hierarchy condition is no longer satisfied , so t hat P^ depend 



also on #12. For these cases we have performed a x squared fit of the three up-down ratios ( |3.5c ) to $23, 6*13 and 



#12, for various values of the mass splittings, and have found acceptable regions of parameter space. Results 
are shown in figure for the case that all Sij = for downward going neutrinos, while all Sij = 0.5 for 
upward going neutrinos and <p = (no CP violation). An equivalent fit would be obtained for cf> = ir and 
#23 — * 7T — #23- The (relatively small) asymmetry of fig.s ^b under #23 — > tt — O23 shows the dependence on <f> of 
the SuperKamiokande data considered here. 

A comparison of figure with figures [| and ^| shows under what conditions this large Amf 2 scheme gives 
consistency. If all solar measurement techniques are correct, then, from figure [| 6> 13 is small and 812 = 10°-h20°. 
Figure ||b then shows that #23 is centred on 45° ± 25°, the range around #23 = 135° being equivalent for any 
4> since 6* 13 is small. Figure || shows that solar model inspired fits to data from two solar techniques at large 
Am 2 2 allow larger ranges of 6*12 and 6*13, and these become correlated with #23 via figure rab. 



5 Large — > v T Mixing For Atmospheric Neutrinos 

The pattern of masses and mixings suggested by the previous considerations show peculiar features, especially 
if both the atmospheric and solar neutrino anomalies are accounted for in the minimal scheme of section [3| The 
mass differences are hierarchical. However a large mixing (#23 ~ 45°) is required between the states with the 
largest mass difference. The mixing angle 612 between the states with the smallest mass splitting may be large 
or small. Finally, if Am 2 tm > 2 ■ 10~ 3 cV 2 , i.e. in the CHOOZ range, the third mixing angle must be small, 
#13 < 13°. Therefore it looks interesting to see which mass matrix could produce this pattern and which flavour 
symmetries can justify it. 



5.1 2x2 Matrix Forms 

As stressed in the introduction, an important consequence of the data on atmospheric neutrino fluxes is the 
need for large mixing angles. Here we study four possible forms of the 2x2 Majorana mass matrix for and 



v r which have a large mixing angle. In subsection 5.2 we study whether these forms can be incorporated in 



3x3 mixing schemes which also give solar neutrino oscillations, and whether 3x3 cases exist which cannot be 
reduced to a 2 x 2 form. In section ^ we study whether these forms may be obtained from flavour symmetries 
of abelian type. 

In a basis with a diagonal charged lepton mass matrix, the Majorana neutrino mass matrix is 



v 2 (C B 



M \ B A 

This is brought into real, diagonal form by the unitary matrix 



(5.1) 



V = R(0)(l e ij (5.2) 

where tan 29 = 2B/(A — C), and the phase a does not affect oscillations. The mass difference relevant for 
oscillations is Am 2 = (A + C)y / (A — C) 2 + 4B' 2 . The coefficient v 2 /M is motivated by the see-saw mechanism, 
with v the electroweak vacuum expectation value and M the mass of a heavy right-handed neutrino. 

There are four possible forms of this matrix which give and these are shown in Table 1. In cases (1) 

and (2) the entries are all of order unity; in the generic case they are unrelated, while in case (2) they are related 
in such a way that the determinant is suppressed. We discuss how such a suppression can occur naturally via 
the seesaw mechanism in the next section. Case (3) has one of the diagonal entries suppressed, which, however, 
does not follow from a simple symmetry argument. For cases (l-j-3), taking Am 2 = 10 _3 eV 2 , one finds 

M = (1 -T- 3) x 10 15 GeV, (5.3) 
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small entries 


small parameters order unity parameters Am 2 / ( jj) 2 


(1) Generic 

(2) Determinant small 

(3) One diagonal small 

(4) Pseudo-Dirac 


none 

none 
one diagonal 
both diagonal 


none A,B,C ps 1 
none A, B, C = B 2 /A + e »1 
C ps e A, B ps 1 
A, C ps e £ we 



Table 1: TTie /our possible 2x2 matrix forms which give a large mixing angle. 



close to the scale of gauge coupling unification in supersymmetric theories. 

Finally, case (4) has both diagonal entries small, making and v T components of a pseudo-Dirac neutrino. 
This follows from an approximate — L T symmetry, and implies that 6 ~ 45°. This agrees well with data: 



combing and ppc of (3.5c) gives 9 = 45° ± 15°. Of the four possible cases with large mixing angle, it is only 
the pseudo-Dirac neutrino which allows to be the astrophysical hot dark matter, in which case one predicts 
9 = 45° to high accuracy. 

From the viewpoint of atmospheric neutrino oscillations alone, the distinction between cases (1) and (2) is 
unimportant. Since case (3) does not follow from simple symmetry arguments, one is left with two main 2x2 
mixing schemes: the generic and pseudo-Dirac cases. 



5.2 3x3 Matrix Forms 

There are many possibilities for 3x3 neutrino mixing giving P w ps 0.5, with oscillation primarily to v T . In 
general two independent frequencies and three Euler angles are involved. 

For the case that the oscillation is dominated by a single frequency, the possibilities may be divided into two 
classes: "2 x 2-like" and "inherently 3 x 3." The 2 x 2-like cases are just the four discussed in subsection 5.1, 
with 6*12,13 small. Even though Am| 3 may not be the largest Am 2 , it is the only one which causes substantial 
depletion of v^. More interesting are the inherently 3x3 cases, for which there is no 2 x 2 reduction. 

Consider the case 



v 

M 



' B 




B 


:) 


K A 






hi = 



(5.4) 



This is diagonalized by V — i?23(6 ) 23)^i2(6'i2 = 45°) giving a Dirac state of v e married to 
The mass eigenvalues are (M, M, 0), which, from the viewpoint of oscilla tions are equivalent 
to (0,0, M). Hence, one immediately sees that the oscillation probabilities are given by ( 3.3c ) with 6*13 



with A, B ps 

C23^V + s 23l / - 

0: 

Pp.T = sin 2 26 , 23S , 23 has th e fo rm of a 2 x 2 oscillation, even though the mass matrix has an inherently 3x3 



form. This arises because (5.4) is governed by the symmetry L e — — L T , which allows 
v e from oscillating. 



but prevents 



We claim that (5.4) is the only inherently 3x3 form for 



v T at a single frequency, as we now show. An 
inherently 3x3 form must have large entries outside the 2x2 block in 23 subspace. The three possibilities arc 
11, 12 and 13 (and their symmetric). None of these entries work alone, even coupled to any structure in the 23 
block: cither one gets two comparable frequencies or one docs not get i/n — > v T . The same is true for 11 + 12 or 



11 + 13, again possibly together with any 23-block. Since 11 + 12 

only case remaining is 12 + 13, with a relatively negligible 23 block, i.e. the 3x3 form in (5.4) 



13 leads to two comparable frequencies, the 

Basic to this 



conclusion is the assumption of no special relations among the different neutrino matrix elements other than 
the symmetry of the matrix itself (for alternatives see |22]). 



6 Models for both Solar and Atmospheric Neutrinos 

In this section we construct models for the minimal scheme for atmospheric and solar neutrino oscillations, 
discussed in section ||. The mass pattern suggested by this scheme has the hierarchy Anig = Ato 2 2 -C Am 2 tm = 
Am 2 3 . We take the form of the lepton mass matrices to be determined by flavour symmetries (FS) and assume 
that all small entries in these matrices are governed by small flavour symmetry breaking (FSB) parameters. 

The low energy effective mass matrix for the three light left-handed neutrinos can be written as the sum 
of two matrices: tull = Tn a tm + m e, where all non-zero entries of m atm are larger than all entries of m®. 
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The form of m atm is such that there is a large mass splitting: Am 2 tm w lO"' 2 ^ 3 -* cV 2 , and a vanishing Am 2 . 
Furthermore, this matrix must give a large depletion of v^, and, as discussed in the last section, this could occur 
if it has certain 2 x 2-likc or inherently 3x3 forms. Of the two 2 x 2-like forms shown in Table 1, only case 
(2) is acceptable: in cases (1) and (3) the two independent Am 2 are comparable, while in case (4) the second 
independent Am 2 is larger than Am 2 tm . Hence, we arrive at the possibility^: 





m£Z = jj C B (6.1) 
M \0 B A ) 

with A, B m 1 and C = B 2 /A. A reason for the vanishing sub-determinant will be given shortly. 

In the previous section we have proved that there is a unique form for m a t m which is inherently 3x3: 

™atm = T7 | B | (6.2) 

with A,B w 1. 

The oscillation angles in the leptonic mixing matrix, V , have contributions from diagonalization of both 
the neutrino mass matrix, 8^, and the charged lepton mass matrix, 6^-\ V(6ij) = V e ' {Oij)V v {6"A. This 
requires discussing also the charged lepton mass matrix. It is not easy to construct an exhaustive list of 
the possible symmetries and their breaking parameters. This is partly because there are both discrete and 
continuous symmetries with many choices for breaking parameters; but is mainly because of a subtlety of 
the seesaw mechanism. Let turu and ttilr be the most general Majorana and Dirac mass matrices of the 
seesaw mechanism allowed by some approximate symmetry. On forming the mass matrix for the light states, 
mii = TnLRTn^ R m/£ R , one discovers that tull need not be the most general matrix allowed by the approximate 
symmetry. This means that one cannot construct an exhaustive list by only studying the symmetry properties 
of ttill — it is necessary to study the full theory containing the right-handed states. 



A casual glance at (S.l) and (3.2) shows that the flavor symmetry we seek, from the viewpoint of AL = 2 
operators, does not distinguish from l T , but does distinguish these from l e . There are many combinations of 
the three lepton numbers L ai and their subgroups, acting on l a , which have this property. As representative of 
this g roup , we c hoo se the combination L e — L M — L T . We find it remarkable that this symmetry group can yield 
both (6.1) and ( |6.2| ), depending on how it is realized. 



6.1 L e — — L T realized in the Low Energy Effective Theory 

In the effective theory at the weak scale, we impose an approximate L e — — L T symmetry, which acts on the 
weak doublets, l e ,fj,,T, and is broken by small FSB parameters, e and e' of charge +2 and -2, respectively, giving 
a neutrino mass matrix: 



„2 



e' 1 1 



™>ll = — ( 1 £ £ | (6.3) 
1 £ £ y 

Hereafter, the various entries of the matrices only indicate the corresponding order of magnitude, allowing for 
an independent parameter for each entry. This texture gives 

v 2 

mi «m 2 «- Am 2 2 ss -^(e + e') (6.4a) 

(6.4b) 





. v2 






m 2 f= 


^ M 


Am 2 2 r 


' M 2 










m 3 s 


- v \ 
* M £ 


Am 2 3 s 


, vi 
Jj M 2 



and 



y 23 



1 we =45°. (6.4c) 



While the texture gives only the order of magnitude of 1 it precisely predicts Q\ 2 to be close to 45°. If the 
FSB parameters e and e' are taken to be extremely small, this becomes an excellent candidate for the case of 
"just so" solar neutrino oscillations, with the prediction that 612 = 45°. However, from figure || it follows that 
this model cannot give matter neutrino oscillations in the sun, which requires sin20i2 < 0.9. There are several 

7 Ansatze of this_type for the neutrino mass matrix, up to small corrections, to describe atmospheric and solar neutrinos are 
contained in ref.s [E3|. 
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contributions to the deviation of sin 20\2 from unity, but they are all too small to reconcile the discrepancy. A 
hierarchy in Am 2 requires e,e' < 0.1, and since sin 2 26^ 2 ~ 1 — (e — e') 2 /8, the deviation of sin 2#5 / 2 from 1 is 
negligible. After performing the B\ 2 rotation, there are small 0(e) rotations in the 13 and 23 planes necessary 
to fully diagonalize tull', these are too small to affect our conclusions. The last hope is that there could be a 
significant contribution to #12 from diagonalization of the charged lepton mass matrix. As mentioned above, 
the diagonalization of the charged lepton mass matrix has to be discussed anyhow. 



Consistently with the symmetry structure of (6.3), the most general form for the charged lepton mass matrix, 



with a structure governed by abclian symmetries is 

m E = \v I £'e £ (6.5) 




when left (right) handed leptons are contracted to the left (right), ei,mE&R- {!,£,,£,') are the relative FSB 
parameters of (tr, hr, e R ) with respect to some other approximate FS, needed to describe the charged lepton 
mass hierarchies, and A is the absolute FSB parameter of TrTl- Here we ignore the fact that non-abclian 
symmetries could modi fy t his form, for example by requiring some entries to vanish. 
Diagonalization of (|6.5| ) leads to 



6*23 ps 1, 9 e 13 ps e', 

Therefore, altogether 



723 



1, 9 12 me + e', 6>i2«45°. (6.6) 



Since sin 26\2 remains corrected only by quadratic terms in e and/or e', we conclude that L e — L M — L T , 
realized as an approximate symmetry of the low energy effective theory, can explain both atmospheric and solar 
neutrino fluxes with a hierarchy of Am 2 , most likely only for the case of "just so" vacuum solar oscillations, 
in which case the scale of new physics, M, is close to the gauge unification scale, and the FSB parameters are 
extremely small: e, e' ps 10 . This result also applies when any approximate FS of the low energy effective 



theory yields (3.2). In view of (3.4), with Amf 3 ps Am 2 tm ps 10~' 2 ^ 3 -' eV 2 , notice that all three neutrinos are 
cosmologically irrelevant. Furthermore, the smallness of the 11 entry of ( |6.3[ ) makes the search for neutrino-less 
2/3-decay uninteresting. 

Comparing the $13 plots of figures I and |, one finds that, with one experiment excluded, the case of 
O12 = 45° becomes allowed for a large range of Am 2 2 , giving another application for this inherently 3x3 form 
of the mass matrix. 



6.2 L e 



L T realized via the Seesaw Mechanism 



The seesaw mechanism |24j allows a simple origin for the vanishing of the 2x2 sub-determinant of (6.1). Consider 
a single right-handed neutrino, N. with Majorana mass M and Dirac mass term vN(cos6v T + svaOv^), where 
9 ml. Integrating out this single heavy state produces a single non-zero eigenvalue in tull — giving ( |6.l| ) with 
A = cos 2 8,B = cos 9 sin 8 and C = sin 2 9, so that AC = B 2 . 

How could this carry over to a theory with three right-handed neutrinos, iV Q ? As long as one of them, N 
with the above mass terms, is much lighter than the others, then it will give the dominant contribution to to^l, 



which will have (6.1) as its leading term. Clearly the key is that there be one right-handed neutrino which is 
lighter than the others, and couples comparably to and v T . 

This can be realized using L e — L^ — L T , with two small FSB parameters e (+2) and e' (—2). The right-handed 
neutrino mass matrix is 

te 1 i\ 

m RR = M 1 e £ (6.7) 
\1 £ e) 

and the Dirac mass matrices of neutrinos and charged leptons are 



m LR 




and niE — Aw | SJe'n £77 77 



(6.8) 



where, in analogy with ( |6.5| ), we have introduced FSB parameters consistent with (6.7) 
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For ease of exposition, let us first consider the case where all the r\ and £ factors are set equal to unity. The 
crucial point is that there is a masslcss right-handed neutrino in the limit e — > 0. Hence, taking e small, and 
doing a rotation in the 23 plane we have 2x2 sub-matrices 



l RR 



1 

M 





1/e 



giving 



M 



rriLR = X'v 



1/e 1/e 
1/e 1/e 



(6.9) 



(6.10) 



where det will = at this order. In a theory with right-handed neutrinos, L e — — L T leads to (6.1) 



Extending the analysis to 3 x 3 matrices is straightforward. The inverse of turu 



1 rr 



e 1 


1 


1 e' 


£ ' 


1 e' 


1 

s 



(6.11) 



shows a pscudo-Dirac structure in the 12 subspace, which is preserved in the light neutrino mass matrix: 

(X'v) 2 




(6.12) 



where we have gone to a basis which diagonalizes the 23 subspace. The parameters relevant for neutrino 
oscillation are 



and 



qe,v 
7 23 



Am 2 



1. 



'13 



23 



1 (X'v) 4 
e 2 M 2 



45°, 



(6.13a) 



(6.13b) 



It is remarkable that L e — — L T has forced a pseudo-Dirac structure in the 12 subspace as in its previous 
realization, again giving #12 near 45°. The crucial difference is that the pseudo-Dirac mass splitting is now a 
higher power in FSB than before 



AW12 

Amj 3 



e 2 (e + e') 



(6.13c) 



rather than e + e' . This allows e and e' to be considerably larger than before, so sin26>i2 < 0.8 is now possible, 
allowing large angle MSW solar neutrino oscillations. In this case the FSB parameters are not very small 
e,e' w 0.3 -j- 0.5, so that the mass of the right-handed neutrinos is still quite close to the gauge coupling 
unification scale. Notice again the cosmological irrelevance of the neutrino masses. For neutrino-less 2/3 decay 
searches (ro u )u ~ e 3 (Am 2 3 ) 1 / 2 < 10~ 2 cV 2 . Finally, e' « 0.1 and A' w 1 can make M exactly coincident with 
the unification scale. 

So far we have only produced models with large 6±2. However L e — — L T realized with the seesaw 



mechanism may also lead to small 612, using the FSB suppression factors in ( |6.8| ). Taking rf <C e' and 77 
in an appropriate 23 basis gives 

(X'v) 2 



M 




so that eq.s (6.13b) and (3.13c) remain valid but 

1, 1 



qe,u 
'23 



13 



and, most importantly 
which can make #12 small. 
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Figure 7: Different combinations of ranges for Am| 3 (dark gray) and Am 2 2 (light gray) discussed in the text. 



7 Conclusions 

The solar and atmospheric neutrino anomalies, strengthened by the recent SupcrKamiokande observations, can 
be interpreted as due to oscillations of the three known neutrinos. However there is still considerable allowed 
ranges of masses and mixing angles that can account for all these anomalies, especially if a cautious attitude 
is taken with regard to the theoretical analysis and/or the (difficult) experiments relevant to solar neutrinos. 
A further major element of uncertainty is related to the relatively large range of values for the mass splitting 
that can account for the atmospheric neutrino anomaly. We summarize our conclusions by considering a set 
of alternative hypotheses, related to these dominant uncertainties, with an eye to the experimental program 
that may lead to their resolution and eventually to the determination of the full set of neutrino oscillation 
parameters. 

A critical value for Am| 3 is around 2 • 10~ 3 eV 2 mainly because for larger values CHOOZ sets a considerable 
constraint on the mixing pattern, but also because (l-j-2) ■ 10 -3 eV 2 is frequently discussed as a typical sensitivity 
limit for various Long-Basc-Linc (LBL) neutrino experiments, like the one from KEK to SK, or the v T appearance 
experiments with a high energy beam from CERN to Gran Sasso or from Fermilab to Soudan. On the other 
end, a value of Amf 2 < 2 • 10 _4 eV 2 , as certainly required by a standard Solar Neutrino Analysis (SNA), would 
make the corresponding oscillation frequency irrelevant to the SK experiment on atmospheric neutrinos. On this 
basis we consider the following four possibilities, none of which, we believe, can be firmly excluded at present. 
They are graphically represented in fig. [7| 

1. Am 23 > 2 • 10 _3 eV 2 and Am 2 2 < 2 • 10 -4 eV 2 . Here a minimal scheme to describe both solar and 
atmospheric neutrinos is required, as discussed in section [§[ with Ato 2 3 ^ > Ato 2 2 . Since Am 2 2 is too 



small to affect atmospheric and/or LBL experiments, in both cases eqs. 3J3 apply. The fit relevant to SK 
is given in fig. ^a, with the further constraint, from CHOOZ, that #13 is small, 6*13 < 13°, and therefore 
$23 = 45°±15°. In turn #i 2 , together with Am 2 2 , will have to be determined by solar neutrino experiments. 
In this alternative, the neatest confirmation of the SK result would come from a v r appearance LBL 
experiment. At the same time, a dominant — * v T oscillation should also lead to a signal in the KEK to 
SK disappearance experiment, with no appreciable v e appearance signal. 

Am| 3 < 2 • l(T 3 cV 2 and Am 2 2 < 2 • 1CT 4 cV 2 . The main difference with respect to the previous case is 
that now #13 is not constrained by CHOOZ, and therefore, from fig. ||a, it can be as large as 45°. This 
implies, from eqs. |3.3| , that the results of both atmospheric and LBL experiments, with low enough 
energies to permit exploration of Am 2 lower than 2 • 10 -3 eV 2 , may be affected by a significant P„ e 7^ 0. By 
the same token, an experiment with low energy v e extending the sensitivity of CHOOZ (e.g. Kam-LAND) 
may show a large signal if $13 is indeed large. In any event P /jr will be significant. Finally, as in case 1., 
decoupling of solar and atmospheric neutrino oscillations implies that #12 can only be determined by solar 
neutrino experiments, with an analysis complicated by $13 being potentially unconstrained (see fig.s ||, 
upper row) 

Ato 23 > 2 • 10~ 3 eV 2 and Am 2 2 > 2 • 10~ 4 cV 2 . This case is possible only if SSM constraints are relaxed 
(fig. |^, lower row) and/or if one of the experimental techniques for solar neutrinos is problematic (fig. |B|). 
However, as discussed in section |I| Am 2 2 must be lower than 2 • 10~ 3 eV 2 , below the CHOOZ range. 
Since, on the other hand, Am 2 tm = Am 23 is in the CHOOZ range, #13 is small and eq.s 4.3 are relevant 
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for atmospheric and LBL experiments. The fit of the present SK results gives #23 = 45° ± 25° (the range 
at 023 ~ 135° being equivalent since #13 is small). Therefore the main difference with respect to case 1. is 
the possibility of a contribution in eq. ( |4.3| ) . While v T appearance in LBL experiments must still give 
a positive signal, P Me could significantly deviate from zero at low enough oscillation frequencies (relevant 
to lower energy LBL experiments or to reactor experiments such as Kam-LAND). The finding of such 
an effect, together with a positive v T appearance signal, would prove, in the three neutrino oscillation 
picture, the inadequacy of the NSA as it is done now. 

4. Am,23 < 2 • 10~ 3 eV 2 and Am 2 2 > 2 • 10~ 4 eV 2 . This is the relatively less constrained case (and also the 
relatively less likely). Here both neutrino squared mass differences are outside of the CHOOZ range, so 
that 613 is unconstrained. Appropriate values of the mixing angles can fit the SuperKamiokande up/down 
ratios of atmospheric neutrinos, as shown in fig. |6^. In this case, the two comparable Am 2 might lead to 
sizeable CP-violating effects if all the three mixing angles are large. 

Measurements by SNO and Borexino will increase the number of independent observational signals of the solar 
fluxes, Si, from 3 to 5; so that, from (2.12) with $cno/ < I' 7 Bc = 0.22, Am 2 2 , 612, 813, | I > 7 Be and can all b G 



determined. This will provide a crucial consistency check between the experimental techniques and the solar 
models. If $13 is found to be large, Am| 3 < 2 x 10~ 3 eV 2 , giving a signal at Kam-LAND, but making it harder 
for LBL experiments. 

In the minimal scheme, with a hierarchy amongst the Am 2 , several years of data from Supcr-Kamiokandc 
will allow a fit to Am 23 ,# 2 3 and #13. Combining with fits to the solar flux measurements, and to LBL and 
Kam-LAND experiments, could allow the emergence of a consistent picture for the two oscillation frequencies 
and the three leptonic mixing angles. 

The variety of possibilities discussed above makes it uncertain which is the relevant neutrino mass matrix 
and, a fortiori, which are the flavour symmetries that might be responsible for it. Nevertheless, focusing on the 
minimal scheme for both solar and atmospheric neutrinos, the peculiar pattern of masses and mixings renders 
meaningful the search for an appropriate mass matrix. As discuss ed in secti on js| on general grounds, two forms 



of mass matrices emerge as being able to describe the data, eq.s (3.1) and (63). Since in the minimal scheme 
Am 2 2 <C 2 • 10 _3 eV 2 , these forms imply that neutrino masses will not give rise to an observable neutrinoless 
double beta decay signal. The combination L e — L M — L T of the individual lepton numbers may play a role in 
yielding both these forms. A common feature of the resulting solutions is that the heaviest neutrino mass is 
determined by the oscillation length of the atmospheric neutrinos, (Am^ tm ) 1 / 2 . As such, the neutrino masses 
are irrelevant for present cosmology. Again quite in general, an increasing separation between the two Am 2 
requires the angle #13 to become increasingly small. 
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